Heteroepitaxial growth of high-K gate oxides on silicon: insights from first-principles 

calculations on Zr on Si(OOl) 
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Metal deposition of Zr an a Si(001) surface has been studied by state-of-the-art electronic structure 
calculations. The energy per Zr adatom as a function of the coverage shows, that Zr forms silicide 
islands even at low coverages. Adsorbed Zr is thermodynamically unstable against the formation 
of bulk silicide ZrSi2. The observation that the islands consist of structural elements of the bulk 
silicide is an indication that silicide grains will form spontaneously. 
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I. INTRODUCTION 

The scaling of the CMOS transistor has been the driv- 
ing force behind the tremendous increase in microproces- 
sor performance observed during the last decades. While 
the problems of the past were dominated by manufactural 
aspects, one now faces the first fundamental physical lim- 
itations as structures in logical devices approach atomic 
dimensions. 

As a rule of thumb, the thickness of the gate-oxide has 
to be directly proportional to the channel length in a 
MOSFET (metal oxide silicon field effect transistor) de- 
vice. In the course of the ongoing miniaturization, also 
the thickness of this insulating layer is being continually 
reduced. By the year 2007, gate-oxides in a transistor 
will approach a thickness of 1.5 nml], which corresponds 
to about ten typical atomic distances. The quantum me- 
chanical tunneling currents through such a thin oxide 
are intolerable, and cause increased power consumption 
and deteriorated switching characteristics of the transis- 
tor. Replacing Si02-based oxides, nowadays employed 
as gate dielectric, is one of the "key challenges" to the 
semiconductor industrypj, which has to find a solution 
within the next 4 to 5 years. 

Employing high-k (= large dielectric constant) oxides 
would allow for a greater physical thickness while preserv- 
ing the electrical properties. Possible candidates have to 
meet an extensive list of requirements such as suffi- 
ciently large band offsets, thermodynamical stability or 
interface quality. 

The Si-Si02 system meets all these requirements in an 
unparallelized wayQ. First attempts to grow alternative 
oxides on Si(001) did not yield satisfactory results for 
a variety of reasons - often originating in the interface 
region. 

In the course of changing to a new gate material, one 
also considers to now grow crystalline oxide layers. It 
has been a paradigm that the gate oxide must be amor- 
phous in order to avoid dislocations and grain bound- 
aries, which provide natural pathways for leakage cur- 
rents and atom diffusion. However, in the very small 



devices used in the near future, the probability for such 
a defect may be sufficiently small. On the other hand, 
epitaxial oxides bear the promise of very low defect con- 
centrations at the interface to silicon, which is crucial for 
efficient device operation. 

The first account of heteroepitaxial growth of a tran- 
sition metal oxide has been presented by a seminal 
work of McKeeQ who succeeded in growing SrTi03 on 
Si(001). For the divalent elements detailed studies have 
been performed 0, 0, 0, H, EJ EJ- As similar under- 
standing does not yet exist for the technological relevant 
early transition metals such as Zirconium. Theoretical 
studies on Zr are limited bulk interface calculations of 
silicates 111). 

The deposition and formation of an interface between 
one of the major contenders for high-k oxides, namely 
Zirconia (Zr02), and silicon has not yet been investi- 
gated with ab-initio simulations. As a first step we per- 
formed state-of-the-art electronic structure calculations 
of the deposition of Zr on Si(001) up to a coverage of two 
monolayers (ML). 



II. COMPUTATIONAL ASPECTS 

We performed first-principles calculations within the 
framework of density functional theory^, 0] using the 
gradient corrected density functional of Perdew, Burke 
and Enzerhof [14| . The electronic structure problem 
has been solved with the projector augmented waves 
method[15j, which uses augmented plane waves to de- 
scribe the full wave-functions and densities without shape 
approximation. The PAW method as implemented in the 
CP-PAW code employs the Car-Parrinello approach [Trj| 
to minimize the total energy functional. 

The core electrons are described within the frozen core 
approximation with the semi-core (4s and Ap) shells of 
the Zr atoms treated as valence electrons. Plane wave 
cutoffs of 30 and 60 Ry for the wave functions and the 
density have been used. 

For metallic system we minimized the Mermin func- 



2 



tional with respect to the occupation numbers which 
yields the Fermi-Dirac distribution for the electrons. An 
electronic temperature of 1000 K was used. The zero- 
Kelvin result has been extrapolated using the method 
suggested by [l7| . 

The reconstructed Si(001) surface is the template for 
our growth studies. It is modeled by a slab consisting 
of 5 silicon layers, where the bottom layer is saturated 
with two hydrogen atoms per silicon. The position of 
the hydrogen atoms and the lowest silicon layer has been 
frozen. Supercells with 16 atoms per layer have been 
used for the calculations of Zr in the dilute limit, 8 atoms 
per layer for all higher coverages. In the first case, the 
distance between two periodic images of an adatom is 
15.36 A. For all surface calculations the k-mesh density 
corresponds to 64 lateral k-points per p(l x 1) surface 
unit cell. 



III. DILUTE LIMIT 

Fig. ^ shows the structure of the reconstructed 
Si(001) surface. Excellent accounts covering the full com- 
plexity of the surface reconstruction can be found in 
literature[HITgj. The resulting structure can, however, 
be explained by the following considerations: While ev- 
ery bulk silicon is 4-fold coordinated, the surface atoms 
lack their upper bonding partners which leaves 2 dan- 
gling bonds per atom, each occupied by one electron. In 
a first step, two neighboring silicon atoms form a dimer 
bond - initially parallel to the surface - which leaves one 
dangling bond per atom. Now both electrons are trans- 
ferred into one dangling bond to reduce the number of 
unpaired electrons to zero. This results in the buckling as 
the silicon with the filled dangling bond prefers a tetra- 
hedral bonding arrangement, while the other prefers a 
planar sp 2 configuration. From the electrostatic point of 
view, the alternating buckling behavior is favorable. Our 
calculations readily reproduce previously reported results 
such as the c(4 x 2) reconstruction and the difference in 
z-coordinate of two silicon atoms within a dimer 1201. 




FIG. 1: Si(001) surface with c(4 x 2) reconstruction. The 
figure represents one supercell of the slab as used as a template 
in our calculations. 

The energy surface of a Zr adatom on Si(001), has been 



obtained by freezing the lateral position of the adsorbed 
atom relative to the slab backplane, while all other de- 
grees of freedom were fully relaxed. We used a grid of 
twelve Zr positions in the irreducible zone of the recon- 
structed p(2 x 1) silicon surface. 

The resulting total energy surface for an isolated Zr on 
top of Si(001) is shown in Fig. [21 Two nearly degenerate 
positions can be identified. One is located in the valley 
right in the middle of two dimers of neighboring rows 
and the other on top of a dimer row between to adjacent 
dimers. A third, local minimum, is located in the valley 
and has an energy 0.30 eV higher relative to the most fa- 
vorable positions. The diffusion is quasi one-dimensional 
with barriers of 0.70 eV for diffusion parallel to the dimer 
rows and 1.63 eV from the valley to the row. 




FIG. 2: The total energy surface of an isolated Zr adatom as 
a function of the lateral position on the Si(001) surface. The 
dimer silicon atoms are located at the bar-shaped maxima on 
this surface. The valleys in the energy map correspond to the 
valleys between the dimer rows. 



IV. COVERAGE OF 0.25 AND 0.5 
MONOLAYERS 

In order to investigate the formation of a continuous 
film, we increased the coverage. The structures with cov- 
erages of 0.25 ML and 0.5 ML have similar energies per 
Zr atom as the dilute limit (see Fig.[3J). A wealth of com- 
plex structures has been found. Here, however, we only 
summarize the main trends of the chemical binding. 

To first approximation Zr prefers a formal 4+ charge 
state - the projected density of states of the Zr d-states 
is located well above the Fermi level. It is well known 
that for transition metal cations the s-and p-electrons are 
located above the <i-states, and have only a minor effect 
on the occupied states. The silicon dimers at the surface 
accept up to two electrons in their dangling bonds. A 
clear structural indication that the dangling bonds are 
filled is the disappearance of the dimer buckling. 
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FIG. 3: The energy per Zr adatom as a function of the 
coverage. Energies are given relative to bulk ZrSi2 and bulk 
Si. 



Further electrons that are supplied at increased cov- 
erage of Zr, occupy the antibonding states of the dimcr 
bond. As a consequence the silicon dimers break up. This 
happens at a coverage of half a ML and above. It should 
be noted, however, that metallic states accept some of 
the electrons so that the number of broken dimers does 
not directly correlate with the number of Zr adatoms in 
the ratio one to one. As a result of the interplay be- 
tween breaking dimers and metallic states we find fairly 
complex reconstructions of the surface structure for in- 
termediate coverages. 



V. SILICIDE FORMATION 

Upon increase of the coverage to a full ML, the energy 
drops by 0.64 eV per adatom (see Fig. |3J). The stable 
structure with ML coverage is shown in Fig.Q] All dimer 
bonds of the surface layer are broken, and the Zr atoms 
occupy the centers of the square array of surface silicon 
atoms in the resulting p(l x 1) reconstruction. Such a 
layer is one structural element of bulk ZrSi2 . 




FIG. 4: The structure of a Zr monolayer. Below the silicide 
layer, subsurface dimers have formed. 

The energy gain is, however, not due to the surface 



geometry but can be attributed to a dimer reconstruction 
of the silicon subsurface. A metastable state without this 
reconstruction is even higher in energy than adsorbed Zr 
at lower coverages. 

During the reconstruction, the silicon atoms in the 
layer underneath the ZrSi surface layer form dimers anal- 
ogous to the bare silicon surface. In contrast to the dimcr 
row reconstruction of the silicon surface, these subsurface 
dimers are not buckled and are arranged in a checker- 
board instead of a row pattern. 

At a coverage of two MLs a second ZrSi layer is formed. 
This configuration is shown in Fig. [S] The second layer 
is nearly identical to the first, but is shifted laterally. 
The double layer is again a structural clement of bulk 
ZrSi2, the structure of which is shown in Fig. EI The ZrSi 
double layer consists of one-dimensional Si zig-zag chains 
separated by Zr atoms. The Zr atoms lie approximately 
in the plane of the upper and lower atoms of the Si chain. 
As for ML coverage the first silicon layer underneath the 
ZrSi double layer exhibits a dimer reconstruction. 




FIG. 5: Structure obtained at a structure of two monolayers 
of Zr. The ZrSi double layer constitutes a structural element 
of bulk ZrSi2 (compare to Fig. 6). Note the dimer reconstruc- 
tion of the silicon layer underneath the ZrSi double layer. 

Despite an additional energy gain from the ML cov- 
erage to a coverage of two MLs of 0.25 eV per adatom, 
even this structure is still 0.67 eV higher in energy than 
bulk ZrSi2. While the atomic process has not yet been 
resolved in every detail, our findings give strong indica- 
tions for the nuclcation of the silicide. 



VI. DISCUSSION AND CONCLUSION 

We presented the results of state of the art ab-inito 
electronic structure calculations aiming at understanding 
the deposition of Zr atoms on a Si(001) surface as it is 
the case in an MBE reaction chamber. 

Our results are summarized in Fig. |3 which shows the 
energy per Zr adatom as a function of the coverage. For 
coverages below one monolayer the the energy is nearly 
independent of the coverage. At a coverage of 1 ML we 
observe a sharp drop in energy by 0.64 eV, followed by 
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FIG. 6: Structure of ZrSi2 bulk. It consists of pure silicon 
layers separated by a ZrSi double layers. This double layer 
exhibits Si zig-zag chains separated by Zr atoms. 



further drops in energy for higher coverages. All struc- 



tures are less stable than bulk silicide. 

Our findings suggest that islands with a local coverage 
of 1 ML or higher are formed even at low coverages. The 
islands contain structural elements of bulk ZrSi2, which 
is more stable than any surface structure. Therefore a 
likely scenario is the formation of bulk silicide grains, 
that disrupt the surface morphology, and are detrimen- 
tal for epitaxial growth. Modification of growth condi- 
tions, such as exposing the surface to a oxygen containing 
ambient, may bypass silicide formation during the first 
growth steps. 
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